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ABSTRACT 


The  wave  making  effectiveness  of  a ship  is  an  important  feature  that 
must  be  fully  recognized  in  developing  a uniformly  valid  slender-ship  theory 
for  the  wave  resistance  problem.  The  usual  results  of  ordinary  slender-ship 
theory  are  invalid  near  the  bow  of  a slender  ship  because  the  theory  fails 
to  take  into  account  the  dominant  effects  of  the  divergent  wave  system  in 
this  region.  A new  approach  based  on  the  matched  asymptotic  expansions 
technique  is  used  to  solve  the  bow  problem.  The  general  solution  of  the 
far-field  problem  includes  the  well-known  wave- free  potentials  plus  the 
potential  due  to  a line  distribution  of  wave  sources  that  extends  outside 
the  limits  of  the  ship  to  guarantee  the  proper  description  of  the  far-field 
flow.  The  bow-near-field  problem  and  the  general  form  of  its  solution  are 
derived.  In  order  to  match  the  two  solutions  in  their  overlap  domain  near 
the  bow,  the  far-field  source  function  is  decomposed  into  two  parts  that 
include  a "slowly"  varying  source  function  and  a "rapidly"  varying  source 
function.  The  slowly  varying  sources  can  be  determined  from  the  solution  of 
the  bow-near-field  problem.  However,  the  rapidly  varying  sources,  which  are 
associated  with  the  disturbances  ahead  of  the  bow  and  consequently  needed  to 
obtain  the  solution  of  the  bow-near-field  problem,  cannot  be  determined  unless 
we  match  the  bow  flow  to  the  middle-body  flow.  The  next  step  in  the  analysis 
will  concern  this  problem. 
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1 . INTRODUCTION 


The  present  investigation  concerns  the  steady  forward  motion  of  a ship 
on  the  free  surface  of  an  inviscid  fluid.  The  exact  formulation  of  the 
problem  is  well-known  (see,  for  example,  Wehausen  1973)  but  generally 
intractable,  because  of  the  nonlinearity  of  the  free  surface  conditions. 

No  mathematical ' progress  is  possible  unless  the  problem  can  be  linearized, 
for  which  purpose  additional  simplifying  assumptions  must  be  made. 

Typical  ship  hulls  are  long  and  slender,  which  is  a logical  feature  to 
exploit  in  the  theoretical  investigation  of  such  problems.  The  first  clear 
attempt  towards  a slender  ship  theory  for  the  case  of  steady  forward  motion 
was  made  by  Cummins  (1956)  . Unfortunately , he  did  not  obtain  any  firm  results 
that  could  be  confirmed  experimentally,  and  his  work  was  forgotten  for  many 
years.  Subsequent  investigations  of  the  problem  were  carried  out  independently 
by  Vossers  (1962) , Maruo  (1962) , and  Tuck  (1963) . 

Vossers  reduced  the  linearized  boundary-value  problem  to  an  integral 
equation  for  the  potential  on  the  ship  using  Green's  theorem.  Then  he 
evaluated  the  asymptotic  expansions  of  the  integrals  in  a region  near  to  the 
slender-ship  and  obtained  a first  approximation  for  the  potential  in  this 
region.  Vossers  approach  requires  difficult  approximations  of  multiple 
integrals  that  discouraged  further  developments. 

Tuck  used  the  powerful  technique  of  matched  asymptotic  expansions  in 
his  investigations.  He  matched  the  inner  and  outer  expansions  in  an  "inter- 
mediate" limiting  region  and  obtained  a first  approximation  for  the  potential 
in  this  region.  His  results  are  essentially  similar  to  Vosser's  approximation. 
However,  Tuck  has  made  rigorous  estimates  of  the  error  involved  in  his 
approximations  and  established  the  required  conditions  for  the  validity  of 
the  results.  In  particular,  he  found  that  the  sectional  area  curve  must 
have  an  absolutely  integrable  first  derivative,  a continuous  second  deriva- 
tive and  a piecewise  continuous  third  derivative. 

It  is  rather  obvious  that  ship  hull  forms  do  not  satisfy  these  stringent 
requirements,  and  Tuck  postulated  that  the  conditions  can  be  relaxed  if  the 
results  are  interpreted  as  approximations  rather  than  limiting  values.  In 
particular,  he  assumed  that  discontinuities  in  the  second  derivatives  of 
the  sectional  area  curve  could  be  tolerated.  This  assumption  was  taken  for 
granted  in  subsequent  applications  of  the  theory,  even  though  Tuck  noted  that 
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"the  complexity  of  the  argument  needed  even  to  prove  sufficiency  (of  the 
conditions)  suggests  that  relaxation  of  the  above  conditions  is  unlikely." 

In  the  next  section,  we  will  show  that  the  full  force  of  Tuck's  conditions 
are  needed  in  order  to  justify  the  validty  of  his  results. 

Maruo  investigated  the  wave  resistance  of  a slender  ship  by  evaluating 
the  asymptotic  expansion  of  Michell's  integral  when  the  draft  becomes  small 
and  of  the  same  order  as  the  beam.  He  obtained  the  same  resistance  formula 
as  Vossers  and  Tuck  found,  but  with  some  additional  end-effect  terms. 

Although  Vossers  and  Maruo  did  not  state  any  explicit  conditions  for 
the  validity  of  their  procedures,  the  similarity  between  their  results  and 
Tuck's  results  suggest  the  possibility  that  their  analyses  are  subject  to 
the  same  conditions  obtained  by  Tuck.  In  fact,  the  results  of  these  different 
theories  become  identical  if  the  above  stated  conditions  are  satisfied. 

Unfortunately,  Maruo 's  work  has  caused  a major  set  back  in  the  develop- 
ment of  the  theory,  because  it  led  to  the  belief  that  slender-ship  theory 
results  can  be  considered  as  a special  case  of  thin-ship  theory.  This  con- 
clusion is  based  mainly  on  the  assumption  that  Maruo 's  analysis  is  valid  for 
all  ship  forms.  There  is  no  doubt  that  Maruo 's  results  are  acceptable  if  the 
above  stated  restrictions  on  the  hull  form  are  fulfilled.  But  for  general 
hull  forms,  we  may  use  simple  arguments  to  show  that  Maruo 's  analysis  becomes 
questionable  altogether.  From  purely  mathematical  considerations,  the 
asymptotic  approximation  of  an  integral  when  one  of  its  parameters  becomes  small 
should  differ  from  the  exact  integral  by  an  error  term  that  becomes  small 
when  this  parameter  is  small.  However,  it  turned  out  that  predictions  of  the 
wave  resistance  from  Michell's  integral  differ  appreciably  from  Maruo 's 
asymptotic  formula  for  ordinary  hull  forms  with  small  drafts  (see,  for  example, 
Kotik  and  Thomsen  1963).  In  fact,  the  asymptotic  formula  predicts  negative 
wave  resistance  at  the  low  and  high  speed  limits!  This  simply  means  that 
Maruo 's  results  cannot  possibly  represent  the  asymptotic  behavior  of  Michell 
integral  for  ordinary  hull  forms.  Actually,  the  shortcomings  of  Maruo 's 
analysis  can  be  attributed  to  the  singular  limits  of  the  treated  integrals  as 
was  shown  by  Ursell  (1960) . We  may  conclude  now  that  the  above  statement 
concerning  the  relation  between  slender-ship  theory  and  thin-ship  theory  is 
basically  unfounded  for  ordinary  hull  forms. 


Apparently,  existing  slender-ship  theory  for  the  wave  resistance  does 
not  correspond  to  the  physical  problem  of  interest.  Because  actual  ship 
hull  forms  are  not  "smooth"  enough  to  fulfill  the  stringent  requirements 
of  the  theory,  especially,  near  to  the  ends.  For  this  reason,  investigations 
of  the  end-effects  rendered  necessary  in  order  to  overcome  the  shortcomings 
of  the  theory.  In  this  respect,  Maruo’s  end  effect  terms  are  considered 
incorrect  here,  because  we  believe  that  his  analysis  is  wrong  if  the  ends 
are  not  smooth  in  the  first  place. 

The  first  successful  attempt  towards  a study  of  the  end-effects  was  made 
by  Ogilvie  (1972).  He  assumed  that,  near  the  bow,  gradients  in  the  longitudi- 
nal direction  are  still  small  compared  to  transverse  derivatives , but  not  to 
the  same  extent  as  slender-ship  theory  assumes.  He  considered  the  particular 
case  in  which  the  resulting  boundary-value  problem  for  the  flow  in  the  near 
field  satisfies  the  two  dimensional  Laplace  equation  in  the  transverse  planes, 
together  with  the  "full"  linearized  free  surface  condition  that  is  familiar 
from  classical  thin-ship  theory.  In  order  to  obtain  a unique  solution  for 
the  problem  without  additional  information  from  the  far-field,  Ogilvie 
assumed  furthermore  that  all  disturbances  vanish  ahead  of  the  bow  and  at 
large  side-way  distances  from  the  body.  He  solved  the  problem  for  the 
special  case  of  a fine  wedgelike  bow  by  satisfying  the  body  boundary  condition 
on  its  centerplane.  Daoud  (1975)  extended  the  solution  to  the  case  when  the 
body  boundary  condition  is  satisfied  on  the  exact  body  surface. 

Perhaps  the  most  encouraging  results  concerning  slender-ship  theory  are 
those  obtained  by  Ogilvie  (1972).  He  compared  predictions  from  his  theory 
with  measured  wave  heights  along  the  surface  of  a wedgelike  bow  and  obtained 
satisfactory  agreements.  However,  the  theory  is  valid  only  in  a small  region 
near  the  bow  and  it  ignores  completely  the  interaction  effects  between  the 
flow  in  this  region  and  everywhere  else. 

Thus,  it  is  necessary  to  incorporate  Ogilvie's  bow-flow  problem,  which 
is  basically  a near-field  problem,  into  the  framework  of  a "uniform"  slender- 
ship  theory  for  the  following  two  reasons.  First,  the  theory  will  enable  us 
to  determine  the  correct  far-field  source  distribution  in  the  bow  region  so 
that  we  may  obtain  more  realistic  predictions  of  the  wave  resistance.  The 
second  reason  concerns  the  solution  of  the  bcw-near-f ield  problem,  where  the 
theory  can  provide  an  estimate  on  the  initial  disturbance  at  the  bow  and  the 
appropriate  radiation  condition  in  the  transverse  planes. 


The  present  investigation  concerns  the  development  of  a "uniform"  slender- 
ship  theory  for  the  case  of  steady  forward  motion.  The  term  "uniform” 
implies  here  that  the  theory  takes  into  account  the  influence  of  the  "end 
effects",  especially  at  the  ship's  bow.  However,  it  is  important  to  realize 
that  the  end  effects  associated  with  slender-ship  theory  are  conceptually 
different  from  those  arising  from  slender-bodies  in  unbounded  fluids.  Generally, 
the  end  effects  in  unbounded  fluids  are  due  mainly  to  violation  of  the 
slenderness  assumption,  as  in  the  case  of  a blunt  nose  or  a stagnation 
point  flow  (see  Van  Dyke,  1964).  But  in  slender-ship  theory,  those  effects 
are  more  closely  related  to  the  nature  of  the  wave  system  produced  by  the 
ship  on  the  free  surface.  Details  of  this  behavior  will  be  discussed  in  the 
next  section. 

Apparently,  the  wave  making  effectiveness  of  a slender  ship  is  a 
distinguishing  feature  that  must  be  fully  recognized  in  developing  the  theory. 

For  this  reason,  some  of  the  established  techniques  of  slender-bodies  in 
unbounded  fluids  could  lead  to  the  wrong  conclusions  if  applied  to  ship 
problems.  In  particular,  representing  the  far-field  flow  of  a ship  by  a line 
of  appropriate  singularities  along  its  axis  does  not  guarantee  a true 
description  of  the  wave  system  generated  by  the  actual  ship.  The  general 
form  of  the  far-field  solution  is  discussed  in  Section  3. 
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2.  FAILURE  OF  ORDINARY  SLENDER-SHIP  THEORY 

It  has  long  been  known  that  experimental  confirmation  of  slender-ship 
theory  results  for  the  wave  resistance  are  even  less  encouraging  than  for 
Michell's  integral.  This  has  led  to  the  belief  that,  even  for  ships  of  small 
draft-length  ratio,  the  vertical  distribution  of  sources  represented  by 
Michell's  integral  is  important.  The  failure  of  the  theory  has  therefore 
been  attributed  to  the  fact  that  the  far-field  sources  are  located  right  at 
the  undisturbed  free  surface. 

The  wave  making  effectiveness  of  a line  distribution  of  sources  placed 
at  the  free  surface  was  studied  by  Ogilvie  (1977).  He  considered  the  special 
case  of  a line  distribution  that  gives  the  same  wave  resistance  as  Michell's 
integral.  When  the  analysis  was  applied  to  the  case  of  a strutlike  body, 
the  resulting  source  function  seems  to  smooth  out  the  discontinuous  behavior 
predicted  by  ordinary  slender-ship  theory  at  the  corners  of  the  body.  The 
smoothing  effect  takes  place  over  a relatively  short  region  and  the  resulting 
source  function  extends  outside  the  limits  of  the  body.  Nevertheless,  the 
wave  resistance  which  can  be  computed  from  the  source  distribution  depends 
critically  on  that  small  amount  of  rounding  off  predicted  by  the  analysis. 

Now  one  might  argue  that  the  shortcomings  of  slender-ship  theory  cannot 
be  blamed  on  placing  the  far-field  sources  at  the  free  surface,  but  rather 
on  the  wrong  prediction  of  their  distribution.  However,  this  is  an  over- 
simplification of  the  problem,  since  the  difficulty  is  more  closely  related 
to  the  peculiarities  of  the  flow  caused  by  a line  distribution  which  is 
insufficiently  smooth.  We  note  that  this  problem  is  not  pronounced  as  much 
in  thin-ship  theory,  because  the  sources  are  distributed  over  the  centerplane 
of  the  ship  and  integration  in  the  vertical  direction  alleviates  the  effects 
of  any  discontinuities. 

In  order  to  understand  the  nature  of  the  difficulties  associated  with 
slender-ship  theory,  we  considered  the  flow  caused  by  a line  of  sources  placed 
on  the  free  surface  in  Appendix  I.  The  investigation  of  the  Appendix  concerns 
mainly  the  asymptotic  behavior  of  the  flow  near  to  the  singular  line.  A 
brief  discussion  on  the  results  of  the  Appendix  is  presented  here. 

The  velocity  potential  for  this  problem  was  written  in  a form  which  j 

separates  the  local  disturbance  type  of  behavior  from  the  wave  effects.  The 
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asymptotic  approximation  of  the  local  disturbance  terms  is  independent  on 
the  smoothness  of  the  source  distribution,  but  requires  the  source  function 
to  be  continuous.  On  the  other  hand,  the  behavior  of  the  wave  term  near  to 
the  singular  line  is  more  complicated  and  highly  dependent  on  the  smoothness 
of  the  source  function.  This  term  was  treated  in  such  a way  that  enabled 
us  to  distinguish  between  the  transverse  and  divergent  wave  systems.  It  was 
found  that  if  the  source  function  has  a continuous  first  derivative  the  con- 
tribution of  the  divergent  wave  system  becomes  negligible  near  to  the  singular 
line,  and  the  usual  results  of  ordinary  slender-ship  theory  are  valid,  as 
expected  in  this  case.  However,  if  the  first  derivative  of  the  source 
function  is  discontinuous  at  some  point  Xg  , the  contribution  of  the  divergent 
wave  system  dominates  the  flow  field  in  a small  region  where  (x-xQ  ) =0  ( j y | , 

but  it  becomes  negligible  again  when  (x-Xg)=0(l)  . 

It  is  now  obvious  that  the  full  force  of  Tuck's  conditions  as  discussed 
in  the  introduction  are  needed  in  order  to  ensure  the  validity  of  his  results. 
But  it  is  unlikely  that  actual  ship  hull  forms  can  satisfy  the  stringent 
requirements  of  the  theory,  especially  at  the  ship  bow.  Therefore,  the 
failure  of  the  theory  can  be  attributed  simply  to  the  fact  that  it  does  not 
correspond  to  the  physical  problem  of  interest.  In  particular,  the  theory 
fails  because  it  does  not  take  into  account  the  dominant  effects  of  the 
divergent-wave  system  in  the  bow  region.  For  this  reason,  we  expect  the 
theory  to  give  unreasonable  predictions  of  the  divergent  wave  system,  since 
it  ignores  their  effects  in  the  first  place. 

These  conclusions  may  be  confirmed  by  investigating  the  results  repro- 
duced in  Figure  1 from  Tsutsumi ' s paper  (1978).  The  Figure  shows  a comparison 
between  the  measured  wave  spectrum  of  the  model  with  predictions  from  Michell 
theory  and  slender-ship  theory.  The  sectional  area  curve  of  the  model  has  a 
continuous  first  derivative  everywhere  that  tends  to  zero  at  the  ends , but 
it  can  be  shown  that  the  second  derivative  of  the  curve  has  a square  root 
singularity  at  both  ends.  Thus  the  requirements  of  ordinary  slender-ship 
theory  are  not  fulfilled,  although  its  predictions  are  bounded  because  the 
source  distribution  is  continuous.  According  to  our  previous  discussion,  we 
expect  the  theory  to  give  wrong  predictions  of  the  divergent  wave  system,  and 
that  is  exactly  what  the  Figure  shows. 
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Figure  1.  Wave  spectra  of  M.No.  331  at  Fn  = 0.319 
(from  Tsutsumi  1978) 


Moreover,  the  Figure  indicates  that  predictions  of  the  divergent-wave 
spectrum  from  Micnell  theory  and  slender-ship  theory  have  similar  behaviors 
(their  extremes  occur  at  the  same  points)  which  deviate  appreciably  from  the 
behavior  of  the  experimental  results  at  large  values  of  9 . This  suggests 
the  possibility  that  Michell's  theory  suffers  the  same  ills  of  slender-ship 
theory,  but  to  a lesser  degree  because  of  the  vertical  distribution  of 
sources.  However,  generalization  of  this  conclusion  is  beyond  the  scope  of 
the  present  work. 

According  to  the  results  of  Appendix  I,  the  effects  of  the  divergent-wave 
system  in  the  near-field  are  confined  to  a small  region  of  o(£  /_)  along 
the  ship,  where  s is  the  slenderness  parameter.  In  principle,  this  region 
is  centered  at  points  along  the  hull  where  the  conditions  of  ordinary 
slender-ship  theory  are  not  fulfilled.  However,  these  conditions  are 
generally  violated  at  the  bow  of  actual  ships  and  the  term  "bow  flow"  will  be 
used  here  to  designate  the  nature  of  the  flow  in  such  regions . The  subsequent 
investigation  deals  mainly  with  the  solution  of  the  problem  at  the  ship's 


* 


bow,  but  the  analysis  can  also  be  applied  to  other  points  along  the  hull 
after  some  modifications. 


In  the  bow  region,  the  asymptotic  behavior  of  the  far-field  solution  is 
characterized  by  the  parameter  <x2/|y|  =0(1)  , where  x represents 

distances  measured  along  the  ship  from  the  bow  and  y represents  transverse 
distances.  The  appropriate  inner  expansion  may  be  obtained  in  this  case  by 
approaching  the  body  along  these  characteristic  lines  (see  Cole  1968) . Thus 
in  the  bow-near-field  problem,  where  y = 0(e)  , we  must  require  that 
/<x=0(£1^2)  in  order  to  maintain  the  same  character  of  the  flow  in  this 
region.  Ogilvie  (1972)  used  physical  reasoning  to  obtain  these  requirements 
in  his  study  of  the  bow  flow  problem. 

However,  as  we  proceed  along  the  body  and  x becomes  0(1)  , the 
divergent  waves  will  move  off  to  the  far-field  and  the  transverse  waves 
dominate  the  flow  field  near  the  body.  The  term  "middle-body  flow"  will 
be  used  to  distinguish  this  type  of  behavior  in  which  the  usual  assumptions 
of  ordinary  slender-ship  theory  are  valid. 

It  is  now  obvious  that  a uniform  slender-ship  theory  must  include  the 
above  mentioned  two  near-field  problems.  Each  problem  should  be  matched 
separately  to  the  appropriate  asymptotic  approximation  of  the  far-field. 
However,  to  complete  the  solution  of  the  problem  we  must  match  the  bow  flow 
to  the  middle  body  flow. 
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3.  THE  FAR-FIELD  PROBLEM  AND  ITS  SOLUTION 


We  shall  use  a coordinate  system,  Oxyz,  with  its  origin  fixed  to  the  bow 
of  the  ship.  The  plane  Oxy  coincides  with  the  mean  free  surface,  and  the 
z-axis  is  vertically  upward.  The  ship  is  travelling  with  a constant  speed 
U in  the  negative  x-direction. 

As  usual  in  slender-body  theory,  we  infer  the  existence  of  a small 
parameter  e,  which  is  a measure  of  the  "slenderness"  of  the  ship.  This 
parameter  may  be  simply  the  ratio  of  the  maximum  transverse  dimension  of  the 
ship  to  its  length,  or  it  may  be  the  maximum  slope  with  respect  to  the 
longitudinal  axis  of  planes  tangent  to  the  ship's  surface.  It  is  assumed 
here  that  the  characteristic  length  of  the  transverse  waves  produced  by  the 
ship  is  comparable  to  the  ship's  length  L,  i.e.,  the  parameter  U2/g  = 0(1) 

(we  also  consider  that  L = 0(1)). 

The  complete  outer  expansion  of  the  exact  problem  is  given  by  Tuck  (1963) . 
Accordingly  the  first  order  term  in  the  outer  expansion  of  the  velocity 
potential  satisfies  Laplace's  equation 


1.3x2  + 3y2  + az 

in  the  fluid  domain, 
is 


$ (x,y,z)  = 0 , 

The  boundary  condition  on  the  mean  free  surface 


(3.1) 


z=0 


*xx  + = 0 ' (3-2) 

where  < = g/U2  . The  potential  J must  also  satisfy  the  appropriate  radiation 
condition,  which  is  essentially  a requirement  that  the  waves  lie  behind  the 
disturbance . 

In  the  far-field,  the  ship  shrinks  down  to  a line  and  the  hull  boundary 
condition  disappears  from  the  problem.  Since  there  are  no  other  disturbing 
agents  in  the  fluid,  the  potential  must  be  analytic  everywhere  in  the  lower 
half-plane  except  for  possible  singularities  on  this  line.  For  this  reason, 
it  is  usually  assumed  that  the  far-field  solution  can  be  obtained  by  a 
distribution  of  appropriate  singularities  along  the  axis  of  the  body. 

This  procedure  can  easily  be  justified  for  slender  bodies  in  unbounded 
fluids.  The  general  form  of  the  far-field  solution  in  this  case  can  be 
expressed  in  terms  of  distributions  of  sources  and  multipoles  along  the  x-axis 
(Ogilvie  1977,  p.  100).  These  distributions  are  singular  along  the  x-axis 


■ 


|; 
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and  therefore  cannot  be  permitted  outside  the  limits  of  the  body. 

However,  it  is  not  evident  that  a distribution  of  wave  singularities 
along  the  axis  of  the  ship  could  produce  the  actual  wave  system  in  the  far 
field  of  a slender  ship.  In  fact,  Ogilvie  (1977)  has  shown  that  the  line  of 
sources  which  can  generate  the  far  field  flow  of  a strutlike  body  and  can 
produce  the  same  wave  system  as  predicted  by  thin-ship  theory  must  extend 
along  the  entire  x-axis.  Thus  in  order  to  guarantee  that  the  far-field 
solution  will  present  accurately  the  wave-making  effectiveness  of  the  actual 
ship,  it  is  essential  to  distribute  the  wave  singularities  along  the  entire 
x-axis.  But  this  requirement  contradicts  our  previous  assertion  that  the 
potential  must  be  analytic  everywhere  in  the  fluid  domain  except  on  the 
axis  of  the  ship.  Therefore,  the  far-field  solution  is  determined  according 
to  the  following  rule: 

"we  shall  require  the  solution  of  the  far-field  problem  to  include 
the  most  general  form  of  solutions  to  the  boundary-value  problem 
provided  that  the  asymptotic  expansion  of  these  solutions  near  to  the 
x-axis  is  non-singular  outside  the  axis  of  the  ship." 

Accordingly,  we  shall  construct  the  velocity  potential  in  the  far  field 
by  the  superposition  of  potential  functions  each  satisfying  indentically 
all  the  boundary  conditions  appropriate  to  the  far-field  problem.  Thus  we 
can  write 

$<x,y,Z)  = *(x,y,z)  + JQ  *n(x,y,z>  ' (3-3) 

where  <p  represents  the  potential  due  to  a line  distribution  of  Kelvin 
sources  along  the  x-axis,  and  the  terms  <t>n  represent  wave-free  potentials. 
The  form  of  the  potentials  $n  can  easily  be  obtained  from  the  work  of 
Ursell  (1968) , thus 


00 

<tn  = cos(n0)f  dk  e ikx  a*  <k.)  K_  ( | k j r ) 
1 J n “ 

— OO 


(3.4) 


where  r,9  are  the  polar  coordinates  in  cross-planes  perpendicular  to  the 
x-axis,  such  that,  y = r sin  9 and  z = -r  cos  9 . Here,  Kn  is  the 
modified  Bessel  function  of  the  second  kind  as  defined  by  Abramowitz  and 
Stegun  (1964).  The  functions  a*(k)  represent  the  Fourier  transforms  of 


i| 


i 


L: 


A 
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the  unknown  coefficients  an  (x)  in  the  expansion.  These  functions  must 
satisfy  the  relation 


a* 

2n+l 


+ a* 


2n-l 


= 0 , 


(3.5) 


with  a*^=0  . It  can  easily  be  shown  that  the  wave-free  potentials  satisfy 
Laplace's  equation  (3.1)  and  the  free  surface  condition  (3.2),  while  the 
radiation  condition  is  satisfied  trivially. 

It  is  convenient  to  use  a Fourier  transform  representation  of  the  poten- 
tial for  the  line  of  Kelvin  sources.  One  such  form  is  given  by  Tack  (1963)  as 


' (x,y ,z)  = -2k 


dk  e ilcx  m*  (k)  • 


I exp[-iAy  + zA2+k2] 

lim  | dA  — *.  :v , 

y-K)  J </A2+k2  - (k  + iu/2) 2 


(3.6) 


where  U is  the  so-called  dissipation  factor  of  Rayleigh.  Here,  m*(k) 
represents  the  Fourier  transform  of  the  source  function  m(x)  , 

CO 

ikx  , . 
dx  e m(x) 


m*(k)  = ^ 
2ir 


It  will  be  assumed  that  m(x)  is  continuous  with  piecewise  continuous  first 
derivative  and  absolutely  integrable.  This  assumption  is  needed  for  the 
analysis  of  the  next  section. 

It  is  preferable  at  this  point  to  let  y-K)  in  (3.6)  and  deform  the  path 
of  integration  in  the  cut  A-plane  into  the  proper  indented  contour  C(k)  as 
shown  in  Figures  2a,  2b,  and  2c  according  to  the  value  of  k.  The  poles  of  the 
integrand  are  now  at  the  points  A=±Aq  ' where 

An  = /k2-<2  , for  |k|  > < ; 

A = i /K2-k2  , for  j k | < k 

Ok 

The  expression  for  4>  will  now  be  put  in  a more  useful  form.  For 
y>0  , the  contour  C(k)  may  be  completed  as  indicated  in  Figure  2.  Since 
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the  physical  problem  is  symmetric  with  respect  to  y,  we  need  to  consider 
only  this  case  (y>0)  and  set  y=|y|  in  the  final  results  so  that  they  will  be 
valid  for  all  y.  The  integral  along  the  dotted  semi -circle  can  be  made  as 
small  as  desired  by  taking  the  radius  R sufficiently  large.  Therefore  its 
contribution  vanishes  as  R-*» 

The  potential  4>  can  now  be  expressed  as  the  sum  of  two  parts,  namely 


+ / 


where 


represents  the  outcome  of  the  integration  along  the  two  sides 


of  the  branch  cut  shown  in  Figure  2,  and  $2  is  the  contribution  from  the 
residues.  Specifically,  we  get 


= -4k  dk  e_lkx  m* (k) 


Re  du 


exp[-r(u  sin|8|  + jyu2-k2  cos3)] 


</u2-k2  + i k^ 


k dk  -ikx  „ . 
-4iri  — e m*(k) 

/k2-<2 


exp  -—  (k2  cosy  - ik/k2-<^  sin[d[  ) , (3.8 

< J 

where  the  above  integral  should  be  taken  along  the  upper-side  of  the  branch 


cut  between  k=±<  . 


4.  ASYMPTOTIC  EXPANSION  OF  THE  FAR-FIELD  POTENTIAL  NEAR  THE  BOW 

In  this  section  we  shall  investigate  the  behavior  of  the  potential  $ 
for  points  close  to  the  bow.  As  was  explained  in  section  2,  the  approximation 
that  is  of  interest  here  is  one  in  which  r = Ole)  and  |x|  = Ole1/2)  . 

Let  us  consider  first  the  asymptotic  behavior  of  the  potential  as 

given  by  (3.7).  Substitute  u = ] k ] /X2+r//r  in  the  integral  with  respect 
to  u in  (3.7),  then  change  k to  l//r  afterwards,  hence 


4k 

/r 


di  e-Ux/,/7  m*  (*//?)  Re 


XdX 


q /X2+r  (<X+i | l | ) 


exp 


- 1 £ I (/X2+r  sin  0 + i X cosa) 


(4.1) 


This  expression  is  now  in  a suitable  form  for  the  intended  approximation  in 
which  the  parameter  x//r  can  be  kept  fixed,  i.e.,  x//r  = 0(1)  . At  present, 
the  function  m*(Z//r)  is  assumed  fixed.  Therefore,  we  can  set  r=0  in  the 
integral  over  X and  obtain  the  first  approximation  for  in  the  form 


>1  = - % 1 d*  6 


-i£x//r  „ 


(V^r) 


Re 


dX  e^(.-ili]Xe~||9l)  + E 
kX  + i 2,  1 


(4.2) 


where  denotes  the  error  involved  in  this  approximation.  The  bounds  of 

this  error  term  are  estimated  in  Appendix  II,  part  A,  and  it  shows  that 

E1  = 0(e1/2)  . 

In  the  expression  (4.2),  we  complete  the  contour  of  integration  with 
respect  to  X as  shown  in  Figure  3 and  use  the  residue  theorem  to  evaluate 
the  integral,  thus 

CO  ao 

f « .-(Vf)  • f _S*  cos(|l|Jsin|e|) 


4in<  j -i£<x/£r  , -<£^cos6  7 , , 

— ~ d£  e m*(Ki//r)  e sin (K-i^sin ; 6 j ) , (4.3) 

rr  ‘ 
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where  the  inner  integral  of  the  first  term  is  interpreted  as  a Cauchy 
principal  value  integral. 


Figure  3.  Contour  of  integration  in  the  X-plane 

Since  m(x)  is  assumed  continuous  with  piecewise  continuous  first  derivative, 
we  can  use  the  convolution  theorem  to  evaluate  the  first  term  in  the  above 
expression  (Bochner  1959,  p.  59).  After  the  supstitution  .\=u2/jl|  in  the 
integral  with  respect  to  X we  evaluate  the  1-integration  first,  getting 

» CO 

. 4<  f i , -<u2cosS  , 2 . i,|,  . , i , i , <— 

$ - — : d5  m(C)  j du  e cos  (<u4sin|  6 | ) sin  (ku|  x—,  j /» r)  - 

x /r  J J 

—X>  Q 


4tt<  | , „ -il*x//r 
— — ! dl  e 

/r  J 


m*  (<1/ J?) 


-<12cos9 

e 


sin(<i2sinj 9 J ) . 


(4.4) 


Let  us  consider  next  the  asymptotic  approximation  of  the  potential  :> 
which  is  given  by  equation  (3.3).  After  the  substitution  1 = ^r  k/<  , the 
expression  for  becomes 


. 


| 


-IS- 


4itk:i  f IdSl  -liicx/Sr  „ . „ , /-, 

4>,  ■=-  ~r ==  e m*(<il//r) 

2 /?  J vl2_r 


exp  -k(£2cos9  - i i/$,2-r  sin | 9 |~]  , 


where  the  integral  must  be  taken  along  the  upper  side  of  the  branch  cut 
between  i=i/r 

Again  the  above  expression  can  be  approximated  simply  by  setting  r=0 
while  keeping  x//r  finite.  Thus  we  get 


dl  m*  (<2//r)  sgnU) 

/r 


exp  -k(22cos0  - i £.|  Jl|  sin|  0 j + E2  ' 


where  E2  is  the  error  involved  in  this  approximation.  We  have  also  shown 
in  Appendix  II,  Part  B,  that  E2  = 0(e1//2)  . 

The  approximation  of  the  wave  potential  <J>  in  the  bow  region  can  be 
obtained  by  adding  (4.4)  and  (4.6),  therefore 


4x  f f , -<U2COS0  , 2 • • , 1 » 1 , 

'v  - — d?  m(5)  du  e cos(Ku4sm0)  sin  (ku  | x-ti  | /» r ) - 


4itKi  f -i£<x//r  . . „ , -k£2cos0  , „ 2 

— — dJl  e m*(<il//r)  sgn(2)  e cos  (xl^sgnS ) 

/r  J 


By  using  the  convolution  theorem  on  the  second  term  in  (4.7),  we  get 


,v  - -7=  d£  m(£)  du  e KU  cos®  cos(Ku2sin0)  sin  (ku|  x-5  | /i/r)  - 
/r 


d£  m(C)  | du  e <u  cos0  cos(<u2sin0)  sin  (ku (x-£)//r)  . 


Although  the  two  terms  in  the  above  expression  can  easily  be  combined 
into  a single  term  with  the  integral  over  5 taken  between  the  limits 

and  C=x  , it  is  more  convenient  to  treat  them  separately  at  this  stage. 
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In  fact,  the  first  term  in  (4.8)  behaves  like  (log  ri  as  r+0  and  therefore 
it  is  singular  along  the  entire  x-axis  (this  behavior  can  be  obtained  by 
using  the  Fourier  transform  representation  of  this  term  as  given  by  (4.3)  and 
expressing  the  principal  value  integral  in  terms  of  the  exponential  integral 
as  defined  by  Abramowitz  and  Stegun  (1965)),  while  the  second  term  is 
well  behaved  along  the  x-axis,  as  can  be  shown  from  (4.7).  Therefore 
equation  (4.8)  cannot  be  used  in  its  present  form  to  approximate  the  potential 
ahead  of  the  bow  unless  we  can  remove  the  singular  behavior  of  its  first  term. 

We  shall  suppose  now  that  the  source  function  m(x)  is  composed  of  two 
parts , namely 

m(x)  = ms(x)  + mR(x)  , (4.9) 

where  ms(x)  is  a "slowly  varying"  function  in  the  sense  that  differentiation 
keeps  orders  of  magnitude  unchanged;  while  mR(x)  is  a "rapidly  varying" 
non  oscillating  function  with  the  understanding  that  differentiation  increases 
orders  of  magnitude  by  e-1/2  , i.e.,  mR(x)  = 0(mR£_1//2).  The  physical 
interpretation  of  this  decomposition  can  be  explained  as  follows:  The 
function  ms(x)  represents  sources  that  are  associated  with  longitudinal 
variation  of  the  ship  form  and  therefore  its  derivative  is  proportional  to 
the  rate  of  change  of  the  hull  form  in  the  x-direction  which  is  0(1).  It  is 
also  expected  that  ms  will  vanish  outside  the  limits  of  the  ship.  On  the 
other  hand,  mR(x)  is  associated  with  disturbances  caused  in  the  bow  region 
and  hence  it  will  be  confined  to  a region  of  OCt1/2). 

The  potential  <j)  can  be  written  now  as  the  sum  of  two  parts,  namely 

* = *s  + *R  ' 

where  4> s 4>R  represent  the  contributions  of  the  sources  ms  and  mR 

respectively  to  the  right  hand  side  of  (4.8).  Thus 

x 00 

8lC  | ,,  I , -<U2COS0  , 2 ■ , , , , , ,r-. 

$ = - — dc,  ms($)  I du  e cos(Ku^sinO)  sin  (<u  (x-?)  /vr)  . 

s /r  J J 

0 0 (4.10) 

In  the  above  expression,  we  have  assumed  that  ms(x)  is  identically  zero 
when  x<0  in  order  to  remove  the  singular  behavior  of  (4.8)  ahead  of  the 
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bow.  Consequently,  the  functions  ms  and  mR  may  be  discontinuous  at 
3t=0  , even  though  their  sum  is  continuous  everywhere. 

Let  us  consider  now  the  potential  4>^  > which  may  be  written  as 


2 Z 


(-) 


n+1 


n=0  (2n+l) i [ r j 


2*1 

n+1 

r j 

dC  |x-5 


i2n+l 


mR(5) 


n -ku  cos  6 . ... 

du  u e cos(ku  sino)  - 


4k 

/r 


d? 


mR(C) 


I*  2 

I du  e <u  cos®  cos(Ku2sin6)  sin (ku (x-£) //r)  , 


0 


(4.11) 


where  the  series  was  obtained  by  expanding  the  sine  function  in  the  first 
term  of  (4.8).  Consider  the  integral  with  respect  to  £ in  the  first  term 
of  the  above  expression;  if  we  write 


Mn(x)  = | d£  | x-5 1 2n+1  m^C) 


then  it  can  be  easily  shown  that 
2n+2 

^2Mn  = 2(2n+l)  ! mR(x>  . 
dx 

We  can  use  (4.12)  to  write  M^x)  in  the  form  of  a Fourier  transform 
representation.  Hence 


(4.12) 


, r . m*(k) 

Mn(x)  = 2(-)n+1  (2n+l)  ! j dk  e'lkX 


2n+2 


(4.13) 


After  the  substitution  k=u/»r  in  (4.13),  we  cam  use  it  to  express  '4.11) 
in  the  form 


r r~  ( \ n 

4 00  -i  nx,'  y t • 1 < 

j>_  = p=-  Z.  (n-1)  ! cos(n6)  du  e mf(u/*r) 

R /j-  n=l  j ^ 


r 


4irici  f -iKux//r  . , . , . -ku2cos0  , 7 . . 

— — du  e m*(<u//r)  sgn(u)  e cos  (icu^sinS)  . 

/r  J R 


(4.14) 


The  last  term  is  here  written  in  the  same  form  as  the  second  term  in  (4.7), 
which  is  equivalent  to  the  last  term  of  (4.11)  . Although  the  first  term 
(the  sum)  of  the  above  expression  is  still  singular  ahead  of  the  bow,  we 
will  show  that  this  behavior  can  be  cancelled  by  the  contribution  of  the 
wave-free  potentials. 

Let  us  consider  the  asymptotic  expansion  of  the  wave-free  potentials  in 
the  bow  region.  From  (3.3)  and  (3.4),  we  can  write  these  potentials  in  the 
form 


= E cos (n9 ) 
u n=0 


dk  e‘lkx  b*(k)  (</2|k| )n  Kn(|k|r)  , 


(4.15) 


in  which  we  replaced  the  unknown  functions  a*  by  (</2  k ) b* . Now  the 

n n 


relation  (3.5)  becomes 

< )2 


b*  + b*  =0 


IWJ  b2n+l  " Din  + °2n-l 

In  the  expression  (4.15),  substitute  Z=/r  k and  then  make  use  of  (4.16) 
and  the  well-known  approximations  of  the  Bessel  K-functions  for  small 
arguments.  He  obtain 


(4.16) 


^ -/r  dZ 


-ilx/Jc  b*  ( %/ /7)  (</2  1 1 1 ) 2 log(|k,|/r) 


f 


+ 2 (2n-2)  ! cos(2n-l)9  dk.  e 

2/7  n=1  J 

— 00 


-ikx//r 


bi  . (2.//r)  • 

2n-l 


2n-l 


1 00 

+ — — E (2n) ! cos2n0 
2/p  n=l 


dk.  e 


-iix/tT” 


K 

l2j 


1 2n 


!b*  + r (</2 | A j ) 2b*  .(Z/Sr)' 

zn-i  zn+i 


-19- 


After  neglecting  the  terms  of  0(r1//2)  in  the  above  expression  it  becomes 

<30 

f di  e'Ux//7b*(i/^) 


^ ~=  2h  (n-1)  ! cosnd 
0 2/r  n=l 


where  the  functions  b*  satisfy  now  the  relation 

n 


(4.17) 


b*  = b* 

2n  2n-l 


(4.18) 


We  note  that  (4.18)  represents  now  the  approximation  of  (4.16)  in  the  bow 
region. 

It  is  evident  at  once  from  (4.14)  and  (4.17)  that  if  we  choose  the 

unknown  functions  b*  such  that 
n 


b*  = -8  m*  , 


(4.19) 


which  satisfy  (4.18)  trivially,  then  the  singular  behavior  of  (4.14)  will 
vanish. 

Finally,  we  obtain  the  following  expansion  for  the  total  potential 
$ in  the  bow  region 


$ * - 

/7 


a£  ms (?) 


du  e-tCU  cosS  cos  (KU2sin9  ) sin  [ku  (x-£) //7] 


0 


4nxi 

,7 


du  e-iKUx/^m*(Ku//T)  sgn(u)  cos  (xu2sin6 ) 

R 


(4.20) 


i 


5.  THE  BOW  NEAR-FIELD  PROBLEM 

• 

As  was  explained  in  Section  2 , the  asymptotic  approximation  of  the 
far-field  flow  near  the  bow  is  basically  characterized  by  the  parameter 
<x2/|y|  , which  is  0(1)  . This  suggests  that  in  the  bow  inner  region,  where 
r=0(e)  we  must  require  that  x = 0 ( e 1 ) , in  order  to  maintain  the  same 

character  of  the  flow. 

Thus,  the  asymptotic  expansion  of  the  exact  boundary-value  problem  in 
an  inner  region  close  to  the  bow  can  be  obtained  formally  by  assuming  that 

If  ' If  = 0(fe"1)  ' 0(fe-l/2>  as  e-0  , (5.1) 

where  f represents  any  ilow  variable  in  the  near  field.  This  statement  means 
simply  that  derivatives  in  the  longitudinal  direction  are  still  small  compared 
to  derivatives  in  the  transverse  direction  but  not  to  the  same  extent  as 
ordinary  slender-ship  theory  requires.  The  solution  to  be  obtained  will  be 
valid  in  principle  only  in  a region  close  to  the  bow,  where  r = 0(e)  and 
x = 0(e1//2)  . 

. The  first  approximation  for  the  potential  <j>  (strictly,  the  first  term 

in  the  asymptotic  expansion)  which  is  based  on  the  above  requirements  was 
obtained  by  Ogilvie  (1977).  It  satisfies  Laplace's  equation 

V+*zz  = ° ' (5‘2) 

in  the  fluid  domain.  The  boundary  condition  on  the  undisturbed  free  surface 
z=0  is 

$xx  + K*z  = 0 ' (5.3) 

which  is  the  same  as  the  condition  satisfied  by  the  far-field  problem.  The 
body  boundary  condition  becomes 

3$/3n  = -U  nx  , (5.4) 

• 

where  n = (n^,n2,nj)  is  a unit  vector  normal  to  the  hull  surface  and  directed 
into  the  body;  M is  a unit  vector,  lying  in  a crossplane,  perpendicular 
to  the  contour  of  the  body  in  that  cross  plane. 

In  order  to  solve  the  problem  uniquely,  Cgilvie  assumed  in  addition  that 
there  is  no  disturbance  present  ahead  of  the  bow  and  that  $ = 0(l/r)  as 


. These  conditions  were  imposed  so  that  the  problem  can  be  solved  without 
the  need  of  matching  to  a far-field  solution.  However,  it  is  appropriate 
to  assume  here  that  any  disturbance  caused  by  the  ship  ahead  of  its  bow  will 
be  confined  to  the  bow  inner  region  in  which  |x|=0(e1//2)  . Therefore,  we 
will  require  the  solution  to  satisfy  the  following  conditions: 


lim  $ 


The  radiation  condition  which  may  be  imposed  on  the  solution  as  r-**® 
can  be  determined  from  the  far-field  approximation  of  the  previous  section. 
However,  we  shall  assume  instead  that  |<j>y|,  |<PX|»  and  |<£xx]  are  all  bounded 
as  it-*®  , and  seek  a general  solution  for  the  problem. 


6.  SOLUTION  OF  THE  BOW  NEAR-FIELD  PROBLEM  AND  ITS  ASYMPTOTIC  EXPANSION 


■ 


The  main  purpose  of  this  section  is  to  obtain  a solution  for  the 
boundary-value  problem  stated  in  the  previous  section.  For  general  hull 
forms,  the  solution  which  satisfies  the  given  boundary  condition  on  the  hull 
surface  can  only  be  determined  by  numerical  methods.  However,  it  is  possible 
to  represent  the  solution  by  a series  expansion  with  unknown  coefficients 
that  can  be  determined  later  from  the  body  boundary  condition.  The  specific 
procedure  for  obtaining  the  coefficients  of  the  series  will  not  be  discussed 
here,  since  we  are  interested  at  the  present  time  in  matching  this  solution 
to  the  far-field  problem.  The  method  which  will  be  used  to  solve  this  problem 
is  based  largely  on  similar  works  by  Ursell  (1949,  1950,  1968) . 

It  is  convenient  to  introduce  the  complex  Z-plane,  which  is  defined  as 

„ . -i9 

2 = y + iz  = -ire 

Consider  the  flow  outside  a circular  cylinder,  |z|  = R , that  encloses  the 
body  and  has  its  axis  along  the  x-axis.  Let 

f(x;Z)  = 4>(x;y,z)  + i ip(x;y,z ) 

be  the  complex  velocity  potential  in  this  region  (where  iJj  is  the  stream 
function  conjugate  to  the  velocity  potential  <J>)  . The  free  surface  condition 
(5.3)  becomes 

Re{fxx  + i<fz}=0,onz  = 0.  (6.1) 


The  radiation  conditions  (5.5)  are  now  in  the  form 

lim  Re{f}  = 0 , 

x-*— 00 

lim  Re{fx)=  0 . (6.2) 

x-*— 00 

Also,  from  the  assumed  radiation  conditions  on  $ as  r-*»  , we  shall  require 
that  f be  regular  and  that  | f 2 1 and  |fxx!  be  bounded  outside  |z|  = R . 

Let  us  Fourier  transform  the  complex  potential  f with  respect  to  x,  putting 


30 

f*(X;Z)  = gr  ' f (x; Z)  e j AX  dx, 

2i  j 
.00 


(6.3) 


LI 


where  we  have  used  "j"  as  the  imaginary  unit  to  distinguish  between  the  two 
complex  planes  used  to  solve  the  problem.  When  indicating  "Re"  and  "Im"  for 
the  real  and  imaginary  parts  of  a function,  we  will  always  suffix  an  i or 
j to  the  symbol  to  show  whether  we  mean  the  real  or  imaginary  part  with 
respect  to  i or  j [see  Ogilvie  and  Tuck  (1969)  for  more  discussion] . The 
free  surface  condition  (6.1)  becomes 


\2 

Imi{f*2  + i — f*}  =0  , on  z = 0 . 

Let  us  now  introduce  the  auxiliary  function 

i2 

G*  ( A ; Z)  = f*z  + i — f*  (X ; Z)  , 

^ < 

which  according  to  (6.4)  satisfies  the  condition 


(6.4) 


(6.5) 


Imi  G* (X ;y+i0)  = 0 . (6.6) 

By  the  Schwarz  reflection  principle,  we  extend  G*  analytically  into  the 
upper  half  of  the  Z-plane 


G*(X;y+iz)  = G*(X;y-iz)  • (6.7) 

From  the  assumptions  about  f it  follows  that  the  extended  function  G* 
is  regular  and  bounded  in  the  whole  Z-plane  outside  J Z j = R . Thus  G* 
can  be  developed  in  a Laurent  series  (see  Wehausen,  1974) 

Ci (X)  C?(X)  , . 

G*  (X  ; Z)  = CQ  (X)  + + ~^2 — + , | 2 J > R . (6.8) 

The  condition  (6.6)  implies  that  all  the  must  be  real  with  respect  to  i. 

We  may  now  integrate  the  differential  equation  in  f*  as  given  by 
(6.5)  and  obtain 


f*(X,-Z) 


-iX2Z/< 

e 


Z 

d£  G*(A;5)  elX2?/K 


The  same  value  of  f*  can  be  obtained  by  integration  along  two  paths  not 
enclosing  the  circle  |zj  = R . However,  integration  along  a closed  path 
surrounding  the  circle  |z|  = R does  not  lead  to  a one-valued  function  due 
to  the  poles  in  (6.8)  . Therefore,  we  introduce  a cut  in  the  complex  Z-plane 


extending  from  iR  to  i°°  along  the  imaginary  axis,  across  which  analytic 
continuation  is  not  permitted.  Since  i}>*=Rei{f*}is  an  even  function  with 
respect  to  y it  is  continuous  across  the  cut,  while  its  conjugate  ip*  may 
be  discontinuous.  Denote  by  Ay*  the  discontinuity  across  the  cut;  since 
G*  is  continuous  across  the  cut,  thus 

\2 

Ay*  - A = 0 , 

Z K 

and 


Ay*  = B (X)  eX2z/K  . 

The  objective  now  is  to  find  a source  function  S which  satisfies  the 
boundary  conditions  of  the  problem  and  for  which  the  Fourier  transform  S* 
has  the  same  discontinuity  along  the  cut  as  f*  . Then  the  function  f*  - S* 
will  be  continuous  across  the  cut,  and  the  cut  may  therefore  be  removed. 

Consider  the  pseudo-impulsive  source  placed  at  the  origin,  as  given  by 
Wehausen  and  Lai tone  (1960,  p 495) 


S = 


x 

r 


-2k 


d£  Q(5) 


00 


0 


dk  -ikZ  , 

e sin  /Kk  (x-<i ) 

/ ick 


(6.9) 


In  which  Q(£)  denotes  the  strength  of  the  source.  The  Fourier  transform 
of  S may  be  written  in  the  form 


S*  = -Q* (X ) 


-ikZ 


k-X2/< 


dk 


0 


(6.10) 


where  the  integral  is  now  taken  along  a path  indented  below  the  pole  at 
k=X2/<  [see  Stoker  (1956)  pl78  for  detailed  discussion] . After  a somewhat 
lengthy  transformation,  it  can  be  shown  that  S*  may  be  continued  in  the 
upper  half  of  the  Z-plane  and  that  it  has  a discontinuity  in  its  imaginary 
part  with  respect  to  i across  the  cut  of  magnitude 


Imi{AS*}  = 2ir  Q*(X) 


(6.11) 
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Consider  the  complex  function 


f*(A;Z)  - S*(A;Z) 

which  is  continuous  across  the  cut  if  we  set  Q* (\)=B(\)/2n  in  (6.11);  the 
cut  therefore  may  be  removed  and  the  function  has  a Laurent  expansion  for 
| Z j >R  . Thus 


f*  - S*  = ?.  q*(A)  (iZ)n  + ",  q*  (A)  (iZ)'n  = F*(A;Z)  + F*(A;Z)  . 

n=0  n n=l  -n  12 

(6.12) 

The  coefficients  q^  and  q*^  are  all  real  with  respect  to  i,  since 
is  even. 

It  is  possible  to  determine  F*(A;Z)  explicitly.  Using  (6.10),  we  can 

write 


. A2  O* (A) 

S*  + i — S*  = »-£-*■  , 

Z ic  Z 


from  (6.12)  and  (6.8),  we  get 


dFi  A 2 

1 + i — F*  = CQ  (A ) . 


dZ 


tc  1 


Thus 


F*(A;Z)  = P*  ( A)  e“lA2z/<  - iCQ(A)K/A2  , (6.13) 

where  P*  is  real  with  respect  to  i from  the  symmetry  of  4>*  . 

Since  both  S*  and  F*  satisfy  the  free  surface  condition  separately, 

F*  must  also  satisfy  (6.4)  and  it  becomes 

F2(A;Z)  = mil  q2m(A) f <l2>  + <(2m-l)  (lZ)  } * (6.14) 

Finally,  we  obtain  the  following  expansion  for  f*  : 


f* ( A ; Z ) 


S*  ( A ; Z)  + P* (A ) 


-iA2Z/<  <? 

e + m=l  q*2m 


(A){(iZ)“2m  + 


X2 


< ( 2m- 1 ) 


(iZ) 


-2m+l 


(6.15) 


-26- 


By  taking  the  real  part  with  respect  to  i of  the  above  expression  and  using 
the  inversion  formula  for  Fourier  transforms,  we  get 


I 


<t>  (x;y,z)  = -2k  dS  Q(£) 


dk  kz  r— 

e cosky  sin/<k(x-c,) 

/ick 


00 

+ ! dX  e_1AX  P*(X)  eA“2/<  cos(X2y/K) 
—00 


+ 


q2m(x) 


cos  2m6 

2m 

r 


cr11  (x) 

H2mv  ' cos (2m-l) 9 

< (2m-l)  2m- 1 

r 


(6.16) 


The  first  term  in  the  above  expression  represents  a line  of  pseudo- 
impulsive  sources  placed  on  the  free  surface.  The  limit  of  this  term  as 
y and  z tend  to  zero  is  singular  (see  the  discussion  following  equation 
(4.8)),  but  since  the  potential  must  be  analytic  everywhere  ahead  of  the  bow, 
therefore  Q(x)= 0 when  x<0  . 

For  the  purpose  of  understanding  the  second  term  in  (6.16)  conceptually, 
it  is  convenient  to  make  a substitution  t=x/U,  with  t interpreted  as  a 
time  variable.  Then  the  integrand  represent  a standing  wave  potential  with 
amplitude  proportional  to  P*(X)  and  wave  length  2ttk/X2  . However,  the 
integral  represents  the  potential  due  to  propagation  of  an  initial  disturbance 
P(t)  as  discussed  by  Wehausen  and  Laitone  (1960,  p 507). 

The  sum  term  in  (6.16)  represents  the  so-called  wave-free  potentials. 

These  functions  are  singular  when  r=0  , therefore  we  shall  require  that 

q2m(x)=q2m(x) S°  ' f°r  X<°  ' 

It  remains  now  to  examine  (6.16)  more  closely  and  decide  what  are  the 
unknown  quantities  that  can  be  obtained  from  the  near-field  solution.  First 
of  all,  we  note  that  there  is  no  physical  problem  to  solve  when  x<0  , because 

the  body  ceases  to  exist  in  this  case.  From  the  above  discussions  the  solu- 
tion ahead  of  the  bow  is  determined  by  the  second  term  in  (6.16),  which 
include  the  unknown  function  P(x)  . Also,  it  can  be  shown  that 


■l>(x;0,0)  = P(x)  , for  x<0 


-27- 


This  means  that  the  function  P(x)  is  definitely  related  to  the  disturbance 

caused  ahead  of  the  bow  and  it  can  only  be  determined  from  matching  to  the 

far-field  problem.  Therefore,  the  near-field  problem  can  be  solved  completely 

for  x>0  if  we  know  P(x)  beforehand . Consequently,  the  unknown  quantities 

Q(x)  and  q_  (x)  in  (6.16)  are  entirely  determined  from  the  near-field 
2m 

solution. 

The  asymptotic  approximation  of  the  near-field  problem  as  r-«°  while 
keeping  x//r  fixed  can  easily  be  obtained  from  (6.16).  Substituting 
k=Ku2/r  in  the  inner  integral  of  the  first  term,  and  A = <u/Sr  in  the 
second  term,  then  as  r-*00  we  get 


4k: 

(x;y,z)  ^ ■— 

vr 


A w 

j d£  Q(£)  [ du  e <U  cos®  cos  (icu2sin6 ) sin  (<u  (x~%)//r) 


-iKUX/ /r  . -KU2COSe  , o • 

du  e P*(Ku//r)  e cos(<uzsin6) 


(6.17) 
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7.  MATCHING  THE  BOW  FLOW 


The  inner  and  outer  solutions  are  matched  in  a suitable  overlap  domain 
in  which  x//r  = 0(1)  and  e<<r/L<<l  . The  limit  of  the  outer  solution  in 
the  overlap  domain  is  given  by  equation  (4.20) , and  the  limit  of  the  inner 
solution  in  the  same  domain  is  given  by  equation  (6.17).  Therefore,  the 
matching  can  be  accomplished  by  equating  the  two  expressions. 

Consider  the  first  term  in  both  equations.  Since  the  ^-integration 
has  a variable  upper  limit  x,  therefore  the  first  term  in  (4.20)  is  equal 
to  the  first  term  in  (6.17)  if  we  set 

ms(x)  = j Q(x)  . (7.1) 

It  is  convenient  to  match  the  second  term  in  both  (4.20)  and  (6.17) 
in  the  Fourier  domain.  Thus  we  obtain 

P*(k)  = - 4iri  sgnCk)  m*(k)  . (7.2) 

The  above  relation  between  P*  and  m*  is  equivalent  to  a Hilbert  transform 
relationship  between  P and  mR  in  the  x-domain  (see  Titchmarch  1948,  p. 
120) . Hence 

f mR<?) 

P(x)  = 4 I -? d?  , (7.3) 

J 5-X 

•40 


where  the  integral  must  be  interpreted  as  a principal  value  integral. 

The  matching  of  the  bow  flow  is  now  completed.  However,  we  still  need 
to  find  the  "rapidly"  varying  source  function  itir(x)  in  order  to  determine 
the  function  P(x)  from  (7.3)  and  solve  the  bow-near-field  problem  completely. 
This  will  be  achieved  when  we  consider  the  next  step  in  the  analysis  and 
match  the  bow  flow  to  the  middle-body  flow. 
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APPENDIX  I.  ON  LINE  DISTRIBUTION  OF  KELVIN  SOURCES 


Suppose  the  line  of  sources  lies  on  the  mean  free  surface  z=0,  along 
the  x-axis,  and  that  z decreases  with  depth.  The  y-axis  is  horizontally 
at  right  angles  to  the  x-axis.  The  co-ordinate  system  moves  with  the  sources 
at  a constant  speed  U in  the  negative  x-direction.  The  source  distribution 
is  assumed  to  extend  from  x=-°°  to  x=°°  , and  it  either  vanishes  identically 
outside  a finite  range  or  decay  rapidly  as  | x | -*-00  . We  will  require  the 
source  density  function  m(x)  to  be  continuous,  but  allow  its  first 
derivative  to  be  discontinuous  at  some  points . 


Determination  of  the  velocity  potential. 

The  velocity  potential  for  this  problem  was  obtained  by  Tuck  (1963) , 
but  his  results  are  not  in  a form  suitable  for  our  purpose.  An  alternative 
derivation  is  briefly  outlined  here,  which  reduces  the  results  to  a form 
similar  to  that  of  Peters  (1947) . 

According  to  Tuck  (1963)  , the  solution  of  the  problem  can  be  expressed 
in  the  form 


$ (x,y,z)  = $ - 2 


dk  e 


■ikx 


m* (k)  k2 


dA  exp (-iAy+z/k2+X2) 


/k2+A2 


</k2+A2  - k2 


where 


m*  (k) 


_1_ 

2tt 


f . 

, dx  e 


m (x) 


is  the  Fourier  transform  of  m(x)  , and  K=U2/g  is  the  wave  number  of  a 
plane  wave  of  phase  velocity  U.  The  integrals  in  the  above  expression  must 
be  interpreted  in  such  a way  that  satisfies  the  radiation  condition  as  will 
be  seen  later.  The  term  ^ represents  the  potential  due  to  a line  of 
sources  which  satisfy  the  familiar  rigid  wall  free-surface  condition,  and  is 
given  by 


1 


00 

f m(;) 

J 4 t(x-C)2+r2jI/2 

—CO 


(Al.l) 
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C\' 
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! 


K 


where 


r * (y2+z2)’1/2. 

3y  use  of  the  Fourier  transform  convolution  theorem,  we  can  write  $ as 


* * *1  ~ 7 


d£  m' (£)  G(x-£,  Y»  z)  / 


(A1.2) 


where  m'  (£)  denotes  the  first  derivative  of  the  source  function,  and  3 is 
given  by 


G(x,y,z)  * -4  In 


dX  cos(Xy) 


kdk  _ exp  (-ikx+zi/kZ+X2) 
/k2+X^  </k^+X2  -k2 


When  the  substitution  t=vk2+X2  is  made  in  the  integral  with  respect 
to  k,  it  becomes 


f 


r 


G(x,y,z)  = 4 Im  ; dA  cos(Xy)  | dt 
J J 


exo (-ixvt2-X2+zt) 
t2  - Kt  - X2 


(Al.3) 


Here  the  path  C is  chosen  as  shown  in  Figure  A-l  so  as  to  avoid  the  zeros 
of  t2-<t-X2,  namely 

t1  = j (<  + /k2+4X2)  > X > 0 


t2  = j « - »,,<2+4X2)  < 0 


and  so  as  to  make  the  integral  over  t vanishes  as  x*-“  . 


/ 

/ 

/ 


Figure  A-l.  Indented  path  of  integration  C in  the  cut  t -plane 
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The  contour  of  integration  C in  (A1.3)  may  be  completed  as  indicated 


in  the  Figure  and  the  final  expression  for  g can  easily  be  obtained  by 
following  the  work  of  Peters  (1947).  Thus  without  more  ado,  the  potential 
function  <t>  may  be  expressed  as  the  sum  of  three  parts,  i.e., 

<j>  = ^ + <t>2  + 4>  3 , 

where 


d£  m' (?)  sgn (x-?) 


du  e~U x-^  .Re 


dy 


exp [iu (ycosy+zsiny ) 1 
u + i k siny 


(A1.4) 


and 


4 


x 

* 

d?  m' (?)  Re 

— CO 


CO 

• 

exp{<2  (u2+l)  +iic  [ (x-?)  -uy]  v'u2+l} 

/u2+l 

—00 


(A1.5) 


We  note  that  the  potentials  and  $2  represent  some  kind  of  a local 

disturbance  behavior  that  decays  rapidly  at  large  distances  from  the  sources, 
while  $2  includes  all  the  wave  effects. 


Approximation  of  the  potential  near  to  the  singular  line 

We  shall  evaluate  now  the  asymptotic  behavior  of  each  part  of  the 
potential  <j>  when  both  y and  z become  small.  The  behavior  of  the  term 
;j>  as  r-K>  is  well  known  from  similar  work  on  slender-body  theory  in  an 
unbounded  fluid,  thus 


4 m(x)  log  r - 2 


d?  m' (?)  sgn(x-?)  log  2 | x-? | 


(A1.6) 


The  integrals  involved  in  the  expression  for  the  term  ^ are  bounded 
and  its  first  approximation  near  to  the  singular  line  can  be  obtained  simply 
by  setting  both  y and  z identically  zero  in  (A1.4) , hence 


! 


< 
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<P2  ^ - 2 


dC  m' (£)  sgn (x-C) 


-u  x-C 


du 


/u2  + k2 


= -it  I dC  m'  (5)  sgn(x-£i)  [H»(<jx-?|)  - Y (k  j x-£  | ) ] , 


(A1.7) 


where  Hg (x)  is  Struve  function  and  YQ (x)  is  Bessel  function  of  the  second 


kind . 


The  asymptotic  behavior  of  the  terms  $ and  <j>2  near  the  singular 


line  can  now  be  obtained  by  adding  (A1.6)  and  (A1.7),  therefore 


$1  + <t>2  ^ 4 m(x)  lo9  r “ 2 | d?  m'  (£)  sgn(x-C) 


log  2 | x-£ I + 


H0(K|x-e| ) ~ j yo(kIx_5I > 


(A1.8) 


The  error  involved  in  the  approximation  of  the  above  expression  is  bounded 
uniformly  by  a function  of  r that  tends  to  zero  as  r-*-0  . The  combination 


"Hq-Yq"  occurring  in  (Al.8)  is  well-known  to  be  a monotone  decreasing 


function  of  its  (positive)  argument.  Furthermore,  the  kernel  of  the  integral 
is  bounded,  since  the  "log"  term  precisely  cancels  the  logarithmic  behavior 


of  the  Yg  function  as  C-’-x  . Hence  the  expression  (Al.8)  is  a valid 


approximation  for  the  complete  potentials  the  source  distribution 

is  continuous  and  has  a piecewise  continuous  first  derivative. 

The  approximation  of  the  potential  4>3  near  to  the  singular  line 
depends  generally  on  how  bounded  the  integral  with  respect  to  u in  (A1.5) 
becomes  as  y and  z tend  to  zero.  Ursell  (1960)  investigated  the  behavior 
of  such  integrals  in  his  study  of  the  wave  system  generated  by  a travelling 
pressure  point  in  the  vicinity  of  its  track.  He  showed  that  the  limits  of 
the  integral  are  non-uniform;  in  particular,  the  wave  amplitude  tends  to 
infinity  as  we  approach  the  track  on  the  free  surface  (as  y-*0  on  z=0) , but 
it  is  finite  if  we  approach  the  track  from  below  the  free  surface  (if  z-*-0  on 
y=0) . However,  Ursell's  results  cannot  be  applied  directly  to  our  problem 
because  we  are  dealing  with  a distribution  of  singularities  and  the  results 
will  depend  on  the  smoothness  of  their  distribution. 
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Let  us  consider  the  integral  with  respect  to  u in  the  expression 
(A1.5).  According  to  Ursell's  analysis,  the  asympvitic  behavior  of  this 
integral  is  determined  by  the  form  of  the  oscillating  term  in  its  integrand. 
The  exact  form  of  this  term  in  our  problem  depends  on  the  outcome  of  the 
^-integration  which  involves  the  specific  details  of  the  source  function. 
However,  it  is  possible  to  evaluate  the  significant  form  of  the  oscillating 
term  in  the  integral  over  u by  considering  first  the  asymptotic  behavior  of 
the  integral  with  respect  to  5 for  large  values  of  jui 

Let  us  consider  first  the  case  when  m'  (?)  is  continuous  over  the  entire 
range  of  the  5 -integration,  it  can  be  shown  that 


as 


..  . ...  i<  (x-t)  /u2+l  i m'  (x) 

d<;  m (?)  e pc  + 

</u2+l 


1 


—GO 

jul-*30  . In  this  case  the  integrand  reduces  to 


u2+l 


exp{<z(u~+l)  - i<y  u/u2+l} 


as 


and  the  integral  over  u is  uniformly  convergent  as  y and  z tend  to  zero. 
Therefore,  the  first  approximation  of  becomes 


^ 4 


x 

f d?  m' (5) 


— JO 
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—00 


du 

/u2+l 


cos  (<  (x-£)  t/u2+l) 


= - 4ir 


dC  m' (C) 


Y0(x(x-O) 


(A1.9) 


If  this  result  is  added  to  equation  (A1.8)  we  obtain  an  expression  for  the 
asymptotic  behavior  of  the  total  potential  $ near  to  the  singular  line 
which  is  identical  to  the  usual  results  of  ordinary  slender-ship  theory , as 
expected  in  this  case. 

On  the  other  hand,  if  the  derivative  of  the  source  function  is  dis- 
continuous at  some  point  x„  , then 


dg  m' {%)  e 


i< (x-£) /u2+l 


k/u--1-!  L 


m’ <x) -dm1 (x.)  e 


i:<  (x-Xq)  v'U2 


“I 


+ 0;--: 
'.U-j 
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as  lul-*®  ; where  Am' (Xg)  denotes  the  jump  in  the  derivative  at  the  point 
Xg  . In  this  case  the  oscillating  term  which  determines  the  asymptotic 
behavior  of  the  integral  with  respect  to  u for  small  values  of  y and  z 
must  include  the  exponential  term  of  the  above  expression.  Therefore,  we 
can  write  the  integral  over  u in  (A1.5)  as  follows: 

30 

du  . . 1 — 

I = Re  i exp [<z  (u^+l)  + i<  (x  -0  vu2+l] 

J vu2+I  0 

—CO 

exp{ix [ (x-Xq) -u|y | ] /u2+l  } . (A1.10) 

This  integral  is  similar  to  the  one  considered  by  Ursell  (1960)  and  the 
subsequent  analysis  will  follow  his  procedure. 

Let  us  apply  the  steepest  descent  method  to  the  integral  (A1.10)  using 
the  oscillating  term  which  becomes  dominant  as  |u|-*»  . The  path  of  inte- 
gration can  be  deformed  into  two  paths  of  steepest  descent  through  the  two 
saddle  points,  which  are  the  roots  of  the  equation 


{ [ (x-Xg)  - u(y|  ] /l+u^  } = 0 


i.e.  the  points 

1 (x-xq) 


, -1/2 

(x-x0)2  „ j 

L y2  J 


which  are  real  if  (x-Xg)2>9y2  . The  oaths  of  steepest  descent  C_  and  C+ 
in  the  cut  u-plane  are  given  by  Ursell  and  shown  schematically  in  Figure  A-2. 


Figure  A-2.  Paths  cut  of  steepest  descent  in  the  cut  u-plane 


Let  us  consider  the  behavior  of  the  integrals  along  C_  and  C+  as 
y and  z tend  to  zero,  provided  that  y/(x-xQ)<<l.  It  is  easy  to  see  that 
the  integrand  tends  to  zero  uniformly  as  |u]-“°  along  C_  , thus 


I 'f 

C -<* 


i3ir/4 


iir/4  ^u2+l 


i<  (x-?)  </u2+T 


This  integral  represents  contributions  from  the  transverse  wave  system  and 
it  can  be  expressed  in  terms  of  YQ-Bessel  function,  hence 


Re  -v  - it  Yq  (k  (x-5)  ) 


(Al.ll) 


We  note  that  the  contribution  of  the  integral  along  C_  to  the  total 


potential 


is  equal  to  the  contribution  from  the  part  ;j>3  in  the  case 


of  a source  distribution  with  continuous  first  derivative.  But  here  we  still 
have  to  evaluate  the  behavior  of  the  integral  along  C+  near  to  the  singular 
line  which  represent  contributions  of  the  divergent  wave  system.  It  is 
obvious  that  the  saddle  point  at  u+  tends  to  infinity  as  |y;-K)  ; thus 
the  contour  C+  moves  off  to  infinity.  To  avoid  this  difficulty,  we  introduce 
the  new  variable  v = u|y|/(x-xQ)  , in  terms  of  which 


-iV4 


i3w/4  Ar2+<S2 


exp[<z (v2+62)/62  + i< (xQ-5) /v2+62/o] 


(A1.12) 


exp  ix  r-T — (1-v)/v2+52  , (A1.12) 

1 y I J 

where  6=| y | / (x-Xq) <<1  . The  saddle  point  is  now  at  v+  = (1+/1-862)  , 

which  remains  finite  as  |y|-H3  . 

It  is  evident  that  the  asymptotic  behavior  of  the  above  integral  depends 
on  the  parameter  < (x-xQ) 2/ | y | . If  the  value  of  this  parameter  is  large 
which  corresponds  to  < (x-Xg)=0 (1)  , then  we  can  apply  the  ordinary  method 
of  steepest  descent  (see,  for  example,  Jeffreys  & Jeffreys  1946,  p 503)  to  get 


(,<  (x-xn)  2 


(x—x  ) ^ ^ 

exp  j k — 1 j exp  [<z/ 462  + i<  (x^-5 ) /26  ] 
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This  result  is  identical  to  the  one  obtained  in  the  case  of  a source  distri- 
bution with  continuous  first  derivative.  Therefore,  ordinary  slender-ship 
theory  results  are  still  valid  in  the  case  of  a source  distribution  with 
discontinuous  first  derivative  at  large  distances  from  the  point  of  discon- 
tinuity. However,  if  the  parameter  < (x-xQ) 2/| y |=0 (1)  the  integral  (A1.12) 
can  be  evaluated  when  S-*0  (the  details  are  omitted) , and  the  potential 
<j>2  becomes 

x x 

* r 2 2 

<f>3  = -4  77  dC  m'  (£)  Yq(k(x-C))  + 4/tT  Im  ftj  d£  m (5)  e ’4  erf  [-ifi  (x-£)  ] 

—00  ~ —00 

(A1.14) 

where  erf  denotes  the  error  function  as  defined  by  Abramowitz  and  Stegun 
(1965) , and 

ft  = [</  (-z+i  | yj  )]1'/2  . 

We  note  that  the  first  term  in  (A1.14)  is  now  small  compared  to  the 
second  term,  because  | ft | =0 (r-1//2)  • Thus  the  contribution  of  the  divergent 
wave  system  dominates  the  flow  field  near  the  singular  line  in  this  case. 
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APPENDIX  II.  ERROR  ESTIMATION 


It  will  be  shown  here  that  under  certain  sufficient  conditions  the 
approximations  given  in  section  4 differ  from  the  complete  potentials  by 
an  error  term  which  is  OCe1/2) . The  regularity  conditions  required  in  the 
analysis  are: 

(i)  The  source  function  m(x)  is  continuous  and  has  a piecewise 
continuous  first  derivative  everywhere. 

(ii)  The  function  m(x)  is  absolutely  integrable  over  the  whole 
x-axis. 


A.  Estimation  of  the  error  Ei 

The  error  term  E^  represents  the  difference  between  the  complete 
potential  given  by  (4.1)  and  its  approximation  (4.2).  Thus  we  get 


-Ux//e 


(x;r,9)  = — | d£  e 
vr 


m 


*(2//r)  Re 


dX 


exp (-i [ 2 | Xcos8 ) 

< X + i 1 2 1 


exp  (- 1 2 \ Xsin  | 9 | ) - ^ exp  (- 1 2 | /X2+r  sin  j 0 j ) 

/X2+7 


(A2.1) 


First,  let  us  investigate  the  x-derivative  of  the  function  E^  . From 
(A2.1)  , the  Fourier  transform  of  the  function  D^sSE^Sx  in  terms  of  the 
variable  2// r can  be  written  in  the  form 

. 4k 


DJ (2//r;r,9)  = -i  — j,  m*  («.//F) 


dX  Re  exp(-i|a  Xcos9) 
<X  + i 2 


exp(-| 4 | Xsin | 8 | ) - 


/X2+r 


exp 


(-|j,|/x2+7  sin  I 9 I ) 


By  use  of  the  first  mean-value  theorem  (see,  for  example,  Jeffreys  & Jeffreys, 
1946,  pp  50),  we  get 


D*  < -i  — — 2 m*  ( V/F)  Re  f 1 -X.Q.c°g9  > 

1 - r ^ kXq  + l | 2 | 


■) 
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f r , 1 

dX  exp(-|2|Xsin  M>  - -==  : exp  ( - 1 2 1 /T^+r"  sin  ] 0 | ) 1 , 
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where  Xg  denotes  the  value  of  X at  which  the  real  part  of  the  function 
attains  its  maximum  value.  The  X-integration  can  be  obtained  analytically, 


| D*  | < — |m*  (i/i/r) 
1 r 


I expj-ij  X,JXncos9 ) I I _ -\l\fr  sin] 8 I 


<XQ  + i\Z\ 


sin  6 


< m*(^/>'/r) 


We  can  write 


| D1 1 < 4k  dk  |m* (k) 


Therefore,  it  follows  that 


= 0(1)  , 


(A2.2) 


provided  that  m*(k)  is  absolutely  integrable,  which  is  a trivial  consequence 
of  the  stated  regularity  conditions  on  m(x).  However,  in  the  bow  region 
x=0(e1'/2)  , consequently 


~r^  = 0 (E  e“l/2)  . 

0X  1 


(A2.3) 


From  (A2.2)  and  (A2.3)  we  finally  obtain 


E1  = 0(el/2)  . 


B.  Estimation  of  the  error  E- 


(A2.4) 


The  expression  for  E2  represents  the  difference  between  the  complete 
potential  in  (4.5)  and  its  approximation  in  (4.6).  Thus 


K 

E2  = 4ifi  — I dJl  e 


-iKx//r  -<12cos6 

m*  (icl/Zr)  e 


sgn(S-)  e 


ixZ  II  sin  0 


l i<Zvl2-r  sinlSl 
/l2-r  S 


(A2.5) 


] 

J 


where  the  integral  is  taken  along  the  upper  side  of  the  branch  cut  between 
H=±/r  . The  error  term  E2  can  be  expressed  as  the  sum  of  two  parts, 
namely 


E2  - E3  + 


where  the  functions  E3  and  E^  may  be  reduced  to  the  following  forms: 


Ej  = 8tk  Im 


. -iicxchy  , ^ , 

dy  e chy  m*(<chy) 


— i 9 — 

exp(-tcrch2ye  ) {exp (-ix | y | chye  Y)-i} 

.it/2 

-8ttk  Re 

' 0 


-Kxcosy  . . , 

dy  e cosy  m*(iccosy) 


exp (-<rcos2ye  ^ 1 9 ) {exp  (-iic  | y | cosy  e *Y)-1} 


and 


E . = 8ttk  Im 


-i<xchy  -y  , ,2  -i|0|. 

dy  e m*  (techy)  e exp  (-terch^ye  1 ') 


tt/2 


f -ixxcosy  -iy  7 — i I Q I 

-8ttk  Re  dy  e m*(tecosy)  e exp  (-tcrcos^ye  ')  . (A2.7) 


(A2.6) 


In  deriving  the  above  expressions,  we  have  made  use  of  the  fact  that 
m*(-k)=m*(k)  since  m(x)  is  real. 

Let  us  consider  first  the  expression  for  E^ . By  use  of  the  first  mean- 
value  theorem,  we  can  show  that  the  absolute  value  of  the  functions  in  the  { } 
brackets  are  << j y J . Thus 


I E3  | < 8ir  1 y [ 


dy  chy  |m*(<chy)|  + 


n/2 

dy  cosy  lm*(«cosy) 


Since  m(x)  is  assumed  to  be  continuous  with  piecewise  continuous  first 
derivative,  hence  |k  m*(k)|  is  certainly  bounded  for  all  k by  the  Riemann- 
Lebesgue  lemma,  and 


(A2.8) 


l«“*(k)  I < |7|- 


for  some  finite  positive  number  M.  Putting  chy=  shy  + e and  use  (A2.8) 
in  the  expression  for  |e3|  , we  get 


f 

[E3l  < 8ttk  l y l | dk  m*(k)  + jl  + — J M 


Therefore,  we  may  conclude  that 


E3  = 0 ( | y | ) . 


(A2.9) 


In  order  to  estimate  the  error  term  E.  , we  shall  consider  first  its 

4 

derivative  with  respect  to  x.  From  (A2.9),  we  obtain 
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Using  (A2.8)  , hence 
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< 8rK:  M 
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1 + T 

3x 

2 

It  follows  that 


3x 


= 0(1)  , 


for  the  same  reasons  that  led  to  (A2.3),  thus 
E4  * 0(e1//2)  . 


(A2.10) 


Finally,  the  estimate  of  the  error  term  E2  can  be  obtained  from  (A2.9) 
and  (A2.10),  hence 


E2  * 0(el/2) 


(A2.ll) 
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